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ABSTRACT 

The charged-Higgs contributions to the decay B — > X s ^f are discussed in the minimal supersym- 
metric standard model at large tan/3. These contributions receive two-loop 0(a s tan/3) corrections 
by squark-gluino subloops, which are nondecoupling in the limit of heavy superpartners and pos- 
sibly large. Their leading parts are already known and were evaluated by using an effective two- 
Higgs-doublet Lagrangian. Subleading corrections coming from higher-dimension operators in the 
effective Lagrangian were ignored, although this is not, a priori, justified when m H ± is not much 
smaller than the typical supersymmetric mass AfgusY- Here, we calculate all subleading terms 
of the 0(a s tan/3) corrections up to 0{{rn\ , w^±/Mg USY ) 2 ), as well as all the exact two-loop di- 
agrams with squark-gluino subloops, beyond the effective-Lagrangian approximation. Comments 
are made on the size of these corrections. 
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Abstract 

The charged-Higgs contributions to the decay B — ■» X s j are discussed in the minimal supersym- 
metric standard model at large tan/3. These contributions receive two-loop 0(a s tan/3) corrections by 
squark-gluino subloops, which are nondecoupling in the limit of heavy superpartners and possibly large. 
Their leading parts are already known and were evaluated by using an effective two-Higgs-doublet La- 
grangian. Subleading corrections coming from higher-dimension operators in the effective Lagrangian 
were ignored, although this is not, a priori, justified when m H ± is not much smaller than the typical 
supersymmetric mass Msusy- Here, we calculate all subleading terms of the 0(a s tan/3) corrections 
up to 0((m 2 , m 2 H± /Mg VSY ) 2 ), as well as all the exact two-loop diagrams with squark-gluino subloops, 
beyond the effective-Lagrangian approximation. Comments are made on the size of these corrections. 

1 Introduction 

The inclusive width of the radiative decays of the B mesons, B — > X s ^, is well described by the short-distant 
processes b — > sj and b — > sg, since nonperturbative hadronic corrections are small and well under control. 
The partonic processes have been evaluated up to the next-to-leading order in QCD. See Ref. pQ for a 
review listing the steps that brought to this achievement, after it was noted that QCD plays a particularly 
important role for the decay B — > X s j J2J. The Standard Model (SM) prediction for the branching ratio 
BR(B — > Xs'f) is, up to today [3], 

BR(B -> A s7 )(SM) = (3.54 ±0.49) x 1CT 4 . (1) 

This average includes results obtained in the various papers of Ref. 0] . The comparison of the SM result 
with the world average jSj of the inclusive branching ratio from recent experiments at Belle CLEO 0, 
and BABAR |S] detectors, 

BR(B -> A s7 )(exp.) = (3.34 ± 0.38) x 10" 4 , (2) 

is rather satisfactory. In the SM, the process b — ► sj, as well as b — > sg, occurs through loops with 
W ± and i-quark exchange, i.e. at the same level in perturbation theory at which physics beyond the SM 
may contribute. The agreement between experimental and SM results for BR(£? — > X s j), therefore, is 
already good enough to constrain exotic contributions. Indeed, it is already routinely used to check whether 
particular directions of parameter space of supersymmetric extensions of the SM are viable or not. 

In the minimal supersymmetric standard model (MSSM), new loop contributions to the two radiative 
decays of the B meson, b — ► sg and b — > 57, come 9 from the charged-Higgs boson, charginos, gluino and 



neutralino. Their contributions are often comparable to or even larger than the SM one. The inclusion 
of QCD corrections to these contributions |1U| . however, is far from having reached the level of precision 
already achieved in the SM, at least for generic regions of the supersymmetric parameter space. The latest 
development in this direction has been the observation that the SUSY-QCD corrections by squark-gluino 
subloops can be significant when tan/3 is very large [111 112) . 

In this talk, we focus on the contribution of the charged-Higgs boson ± and analyze these two-loop 
SUSY-QCD corrections. So far, only the nondecoupling parts of these corrections have been evaluated, by 
using an effective two-Higgs-doublet (2HD) Lagrangian where squarks and gluino are integrated out [TTllr^] . 
Subleading parts of these corrections, i.e. of order ((m^ , m^ ± )/M| USY ) n , generated by higher-dimensional 
operators in the effective Lagrangian, were omitted in these studies, although they could give potentially 
large contributions when m H ± is nonnegligible with respect to the typical supersymmetric mass MgusY- 

Here, we calculate some of these subleading corrections, i.e. up to (9((mf , rn 2 H± /M| USY ) 2 )- We also 
evaluate exactly all two-loop diagrams correcting at 0(a s tan/3) the charged-Higgs contribution to b — > sj 
and b — > sg. This calculation, clearly, encompasses the effective-Lagrangian approximation and includes all 
subleading terms ((m 2 , TOj ? ±)/-Mg USY ) n . I n this talk, after discussing the effective-Lagrangian formalism in 
Section[51 the charged-Higgs contribution to b — > sj and b — > sg in Section|31 we show in Section^Jnumerical 
results for the charged-Higgs contributions to the Wilson coefficients Cj and C$ at the matching scale, 
~ Mw- The evaluation of these same coefficients at low scale, ~ m,b, and therefore of the BR(£? — > X s j), 
becomes, at this point, straightforward. We have, however, refrained from showing results for BR(_B — > 
Xgj) and from extracting exclusion plots for the charged-Higgs boson mass, since the exact calculation 
of the 0{a s tan/3) corrections to the the contribution, or at least the calculation of the subleading 
0(m 2 /MJusy) an d C(( TO t /-^susy) 2 ) terms of these corrections, is not yet available. We comment on the 
outcome of our calculation in Section 0] and we conclude in Section [5] 



2 H couplings to quarks 
2.1 Tree-level couplings 

The MSSM has two Higgs doublets H D and Hu, H D = (H D ,H D ) and Hu = (H$,H$), which break the 
SU(2)xU(l) gauge symmetry through the vacuum expectation values (VEVs) of their neutral components. 
The two VEVs are related to the VF-boson mass as Myy = g\v 2 jl = g\{v 2 D + vfj)/2. Their ratio is 
conventionally called tan/3, tan/3 = vjj/vd- They form the following mass eigenstates: two CP-even scalars, 
(h°, H°); one CP-odd pseudoscalar, A ; the two states of a charged-Higgs boson, H^; and the unphysical 
Nambu-Goldstone modes, (G^, G°). The charged scalars (H ± , G ± ) are related to the charged components 
of the gauge eigenfields Hp v as 

/ G ± \_ ( cos/3 -sin/3 \ / H± \ f . 

{ H± J~{ sin/3 cos/3 J\ H § J ' ( j 

At the tree-level, the couplings of the Higgs doublets Hi(i = D,U) to quarks obey the selection rule of 
the 2HD model of Type II. That is to say, (dn)i = (d,s,b)n couple only to Hu, whereas (ur)* = (u,c,t)n 
couple only to Hu, as shown by the interaction Lagrangian in the gauge eigenbasis: 

C = d B Y d q L -H D -u R Y u q L -Hu + (h.c), (4) 

in which the SU(2)-invariant multiplication of doublets was adopted (A ■ B = eijAiBj, with e 12 = —£21 = 1 
and en = £22 = 0). This constraint is a consequence of supersymmetry. 
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After diagonalization of the Yukawa matrices Y q , Y q — > diag(/i g ), and the breaking of SU(2)xU(l), 
the 6- and i-quark acquire masses 

nib — h b vn — h b vcos(3, m t = h t v\j = htvsm.fi, (5) 
and the H + couplings to the t- and 6-quarks become 

C = V tb h b sin/3 H + i L b R + V tb h t cos/3 H+i R b L + (h.c.) . (6) 

Here, Va is a element of the CKM matrix. Of the two couplings 

g{H + i L b R ) = V tb h b sm(3 = V tb ^- tan/3 

v 

g(H+t R b L ) = V tb h t cos(3 = V tb ^ cot/3 (7) 

v 

the first is greatly enhanced for large tan/3. 
2.2 Couplings up to 0(a s tan/3) 

After supersymmetry breaking, loop corrections depending on soft supersymmetry-breaking parameters 
generate effective couplings of quarks di R to Hjj and of um to Hp, which are forbidden at the tree- level. 
Squark and gluino loops inducing the couplings bnu^H^ and di J t R H^ ) , are shown explicitly in Fig.QJ 

< \ 
~ b R / \ u L h / \ d L 

1 X 1 X ^ 



bn g u l t R g d L 

Figure 1: b R ULH^ and dht R H~D vertices generated by squark-gluino loop, shown in the gauge eigenbasis of 
squarks. 

At momentum scales sufficiently smaller than the typical supersymmetric mass Msusy, these contri- 
butions can be expressed in terms of an effective 2HD Lagrangian (not of Type II anymore) 

C cS = d R Y d qL-{H D -A dR , q Htj) - u R Y u q L -{Hu+A UR>q H c D ) + (h.c), (8) 

where Hp = (—Hp, H^) and Hfj = (— H^* , Hy) , respectively, and where the index q in Ad R , q and A„ Hi9 is 
understood to denote a left-handed quark. The effective couplings Ad R , q and A URtq induced by SUSY-QCD 
loops (as well as Higgsino-squarks loops) are of order 

A dR , q , A UR>q = O = 0(a s Mg VSY ) . (9) 

i.e. do not decouple |131 1141 ITS] in the limit of heavy superpartners, A/susy — * oo. Although sufficiently 
smaller than unity, they can induce large corrections to the couplings to quarks for large tan/3. 

After the breaking of SU(2)xU(l), masses and iJ ± -couplings of the t- and 6-quarks become, at this 
order, 

£(mas S ,mt) = ~h b v cos/3 [1 -\-A bR>b tan/3] bb - h t v sin/3 [1 +A tRit cot/3] tt 

+V tb h b sm(3{l-A bR , t cotl3}H + i L b R + V tb h t cos/3 [1 -A tR>b tan/3] H + t R b L +(h.c.) ,(10) 
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where CP-violating phases in supersymmetric parameters were dropped for simplicity. Notice that the A bR ^ q 
and At Hi9 of Eq. JHj arc now split into A bR .b, A brt j and At R , b , At Rt t- The elements of each of these two pairs 
differ by SU(2)xU(l) breaking effects. This splitting is achieved by including in the effective lagrangian of 
Eq. (JSJ) additional higher-dimension operators. See for example b^q^ ■ H^H^Hu and bn(qL ■ Hd)HjjHu- 

The first line of Eq. i|10|) denotes the running 6-quark mass within the SM, that is to say: 



m 6 (SM) = h b v cos/3 [1 + A bR:b tan/3] . (11) 

The corrections A bRtb tan/3, of 0(a s tan/3) are potentially large ^Sj for tan/3 *5 (A bRtb )~ 1 . As a result, when 
parametrized by m&(SM), the coupling H + ti I bn may significantly deviate |14II15| from the (renormalization 
group improved) tree- level result, 

g (H+izfin) (eff) ~ V tb h bS m0 = V tb mb(SM) J** (12) 

Note that Eq. (|12|) incorporates the resummation of all (a s tan/3)™ corrections to m b 

The corrections to the vertex H + ind,L are factors A-t R ,d tan/3 and also potentially large. The relevant 
loop diagram is shown on the right side of Fig. These corrections give rise to effective couplings 

g (H+i R d L ) (eff) ~ Vtd 1 ^^ (cot/3 -A tR . d ) , (13) 

where the second term in parentheses can very well be of the same order of magnitude of the first. 

For momentum scales nonnegligible with respect to the superpartner masses, the above effective La- 
grangian has to be enlarged to incorporate additional higher-dimensional operators that have momentum 
dependence. For example, the diagram in Fig. [5] generates effective operators of dimension six, such as 



h/ \ dL 

— l x 



ti g d L 

Figure 2: diXhH~D vertices generated by squark-gluino loops, shown in the gauge basis of squarks. 

qLHul^id^qL) ■ Hjj, and in particular, H~j D d,L"l tl {d fl ti,)H'ij , where gets a vacuum expectation value. 
This gives rise to a new set of couplings H+tLj^dL, which in the large tan/3 limit are of size 



g(H+(q)i L (p)rd L (q-p))(eS) ~ V^^-^- A tL , dS m(3 ~ V tdP ^-^^ A th , d . (14) 



, nit (SM) . nit (SM) 

IVI SVSY JW SUSY 

These operators are suppressed by inverse powers of Msusy and decouple in the limit of heavy superpartners. 



3 H + contributions to b — > S7 and b — > sg up to 0(a s tan/3) 

It is well known that b — > S7 and b —> sg are very sensitive to tan/3 |16| . at times even dangerously so |17| . 
However, at the one-loop level, the charged-Higgs contributions to these decays remains fairly independent 
from this parameter, for tan/3 ^ 3. In this regime, the dominant part of the contribution to the decays 
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b — > S7 and 6 — * sg, comes from the diagrams in Fig.|3 where the photon and gluon have still to be attached 
in all possible ways. In this Figure, it is shown explicitly how the factor tan/3 for the vertex tLbRH^ is 
cancelled by the mixing between Hd and Hjj, which brings in a suppression factor sin 2/3 ~ 1/ tan/3. 
Phrased in an equivalent way, the decay amplitude for this diagram becomes insensitive to tan/3, due to the 
cancellation of the factor tan/3 between the tLbnH + and Sz,tnH~ vertices. 



- -X- - v 

HE/' \# c 7 

1 x i 

bn t L t R s L 

Figure 3: b — » sy by charged-Higgs exchange, with helicity flip on the internal fermion line. The Higgs doublet 
constituent of the charged-Higgs state is indicated explicitly. 

However, since the vertex H^s^r is generated at the one-loop level, it is possible to avoid this tan/3 
cancellation, when considering two-loop contributions to b — > sj and b — > sg. The corresponding diagram, 
where again, the photon and the gluon are still to be attached in all possible ways, is shown on the left of 
Fig. El Since also the vertex B^Shth gets generated at the one- loop level, there is another diagram in which 
it is possible to avoid the cancellation of tan/3 of the lowest order contribution. This is shown on the right 
side of Fig.0] Note that, being the vertex H^SLtL of decoupling type, also the corresponding b — > sy and 
b — ► sg contributions decouple in the limit Msusy — > oo. Both classes of contributions coming from the two 
diagrams in Fig. 01 are, nevertheless, of 0(a s tan/3) with respect to the lowest order contribution. They are 
potentially large, without invalidating perturbation theory: their largeness derives from the suppression of 
the lowest order term. 



Br 



Br 



1— X-l 



tL t R g 



SL 



Br 



~^HE 



SL 



Figure 4: Gluino corrections to the H s^tn and H s^l vertices contributing to the decays b sj and 
b — » sg. The photon can be attached to the t, t, sl, and H~ lines, the gluon to the t, t, sl, and g lines. 

Overall, potentially large 0(a s tan/3) corrections to the charged-Higgs contribution to b — > sj (and 
b — > sg) may come from: 

1. mass corrections, 6mj,, to the vertex H + tLbR l|12ll . see Ref. |17|: 

2. proper vertex corrections to the coupling H~SLtR ^2 El as we ll as H~sjJ,l, shown in Fig.0J 

3. corrections to the couplings B~si,tjd and H~ s^t^y, obtained from the two diagrams in Fig.^Jwhen 
the photon is attached to the t- or the s-squarks, and corrections to the couplings H~si,tRg and 
H~s~LtLg, obtained from the two diagrams in Fig. 0] when the gluon is attached to the t- or the 
s-squarks, or the gluino. 

These corrections, induced by squarks and gluino subloops, have been so far calculated by using the 
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effective 2HD Lagrangian of Eq. JSJ), with all superpartners integrated out *. This approach is justified as 
far as all the momenta of the fields in Eq. © are sufficiently smaller than the squarks and gluino. This 
condition is clearly satisfied for the counterterm Snib and its contribution to the H + tLbn vertex, see Eq. (|12[1 . 

For the corrections to the H~s~Ltn vertex, the validity of the effective Lagrangian of Eq. JHJ) requires rrit 
and m H ± to be smaller than Msusy- The condition m H ± <C AisusY, however, is often violated in the case 
of well-known candidates for the mechanism of supersymmetry breaking. Thus, in the resulting models, 
and in general when m#± ^ Msusy, it is possible that the amplitude of the full two- loop diagrams deviates 
significantly from that obtained by making use of the effective 2HD Lagrangian of Eq. @. One may try 
to handle the case of m H ± nonnegligible with respect to Msusy by including in the effective Lagrangian 
higher-order operators, which are suppressed by inverse powers of M| USY . This is equivalent to making an 
expansion of the two- loop integrals corresponding to the diagrams in Fig. 0] with respect to nr^± /M| USY and 
TOt/M| ugY . The operators inducing the vertex H - ^^^ (see discussion at the end of Section |2~2"|) . as well as 
those inducing the vertices H - )}^^, H^SLthl, and H~si,tR.g, H~SLth9, are already of decoupling type, 
i.e. suppressed by inverse powers of M| USY . Nevertheless, additional expansions of the two-loop integrals 
may still be needed. 

A systematic way of making these expansions ad simultaneously keep into account all needed higher- 
dimensional operators is provided by the Heavy Mass Expansion (HME) JHj- Using this technique, we find 
for example that, among the operators of dimension six to be added, it is necessary to include also the very 
same Oy and <D$, 

= T^™b(v)s L <J^b R F^ , 8 ( M ) = ^ mb (^)s L a^T a b R G a ^ , (15) 

which will be part of the effective Hamiltonian used for the calculation of amplitudes of radiative b decays. 
F^ u and G" in these operators are the field strengths of the photon and the gluon, respectively. Details on 
the use of the HME are given in Ref . [20] ■ 

Through the HME technique we have evaluated terms up to 0((m 2 H± ,m 2 

/Msusy) 2 )- Notice that the 

first term of this expansion is, for small values of rriH±, of the order of magnitude of the SU(2)xU(l)- 
breaking corrections to the coefficients of the effective Lagrangian ©; see Ref. [2H - F° r m H ± closer to 
or even larger than Msusy j however, one may question even the validity of this expansion. Of course, 
the truncation up to the 0{(m 2 H± , m 2 /M| USY ) 2 ) can only be justified by comparing with the result of 
the complete two-loop calculation, which clearly goes beyond the effective-Lagrangian approach. We have, 
therefore, also calculated all two-loop diagrams exactly, using methods described in Ref. 22 . In both cases, 
i.e. whether we make use of an expansion or not, we need to calculate the 0{a s tan/3) terms arising from both 
diagrams shown in Fig. Nondecoupling corrections in the effective-Lagrangian approach come from the 
left diagram with photon/gluon emitted only from the i-quark or charged-Higgs boson. All other diagrams 
are decoupling in the Msusy — ► oo limit, and have not been included in previous studies |HH12lE5ll21| . 



4 Numerical results 

We present here results for the charged-Higgs contributions to the Wilson coefficients Ct{^w) and Cs(fiw) 
of the operators O7 and Os at the matching scale ~ [iw, which we have chosen to be M\y- We incorporate 
the corrections to CV(/ivi/) and Cg(pw) discussed in the previous Section. Our normalization of these Wilson 
coefficients is the conventional one, as follows from the definition of the effective Hamiltonian used for the 
calculation of amplitudes of radiative b decays, 

H cK D ~^K V* (CVUWm) + C s (p)0^)) . (16) 

*A similar procedure has been applied to collect higher order corrections of type 0(a tan/3) in Ref. 1181 . 
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Figure 5: Ratios rr,a(jj,w), defined in Eq. $17 f. as functions of m H ±. The supersymmetric spectrum considered 
here is the spectrum I, given in the text. The dotted lines show the goodness of the leading-order approximation of 
the two-loop calculation, the dashed and dot-dashed lines, the goodness of the approximation in which the first and 
the second subleading terms in (m 2 H ±,mt)/Mg VSY are included. 
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Figure 6: Same as in Fig.0 for the spectrum II. 



We show results in which the two-loop diagrams in Fig. 0] are calculated exactly, and results in which 
the off-shell vertices H~si,tn and H~ Shth are treated in an effective-Lagrangian approach, with inclusion of 
leading and subleading terms up to overall suppression factors 1 / (Afg US Y ) 2 ■ We compare these approximated 
results, which we call CV,8(/iwO|approx) with the exact result, denoted by Ct^(i^w)\ exact- To make this 
comparison more transparent, we plot in Figs. |S] and the ratios 

/ s_ Ci (Mw) | approx _ Ci(/lw) |cxact /. _ s f1 -s 

n{fxw) = 7T7 si (* = 7,8), (17) 

^iKP-W J I exact 

as functions of m#± . In these ratios, the mass correction dmb to the vertex H + tLbn cancels out. Therefore, 
we need to specify only a relatively small number of parameters. 

For Fig. we have chosen a heavier superpartner spectrum, called here spectrum I, (m s - t , , m^ 2 ) = 
(700, 500, 450) GeV, the stop-sector mixing angle cos0 t = 0.8, tan/3 = 30, m~ g = 600 GeV, and (i = 550 GeV, 
whereas for Fig. a lighter spectrum is considered: (rris L , , TOj 2 ) = (350,400,320) GeV, cos9 t = 0.8, 
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tan/3 = 30, m,g = 300 GeV, and \i — 450 GeV. This is denoted as spectrum II. As for other input parameters, 
we have used m t (nw) = 176.5 GeV, which corresponds to a pole mass of 175 GeV, and a s (/j,w) = 0.12. 

In both sets of Figures, the dotted lines denote the ratios rf{jj,w) and rs(pw) in which only the leading- 
order approximation is used for the evaluation of the two-loop diagram, i.e. the nondecoupling contribution 
of Ref. JT] ; the dashed (dot-dashed) lines denote the ratios in which the first (second) subleading terms in 
(m(,ra^ ± )/-Mg ugY are also included. 

These Figures show explicitly that for m#± sufficiently smaller than Msusy, the use of the HME, and 
therefore of the effective-Lagrangian approach, is indeed meaningful: the inclusion of an additional term of 
the expansion brings closer and closer to the exact result. For m H ± ~ Msusy, the improvement provided 
by additional terms stops being so evident, at least in the case of the coefficient Cs(/iw)- Going further, it 
becomes quite clear that the results obtained from the HME, and from the effective-Lagrangian formalism, 
cannot be extended to values mjy± ^ Msusy (i.e. strictly speaking outside the range in which they had 
been derived), as one may have hoped. For these values of m H ±, the series of corrections in inverse powers 
of M| USY does not seem to be convergent. 

In the safe region m H ± < Msusy, (that is, safe for the HME), we find that the largest contribution 
to Ct(iiw) comes from the diagram on the left side of Fig. 0] with the photon emitted by the i-quark, for 
Cs(^w) from the two diagrams still on the left side of Fig. with the gluon emitted by the i-quark and 
the gluino. This is true for the exact calculation and for the approximations at different order in 1/M| USY . 
Above this range of values for m H ± , at different orders of approximation, the contributions from different 
diagrams tend to grow differently with ni H ±, producing the crossing points of different lines, visible in 
Fig. Eland partially in Fig. 03 For example, in the approximation in which all terms 0(m 2 ,m 2 H± /Mg VSY ) 
are retained, for a specific value of m_f/±, all terms in l/Mj USY cancel out and the value of the Wilson 
coefficient coincides with that in which only the nondecoupling terms are kept. Similar cancellations of 
terms happen in the other two crossing points. In any case, there is no particular meaning in these points, 
since they appear in a region in which the HME is not well behaved. 

For m H ± sufficiently smaller than Msusy, the approximation in which only the nondecoupling terms of 
the two-loop diagrams are included, is not a bad approximation of the exact calculation, for the coefficient 
Cr(/j,w)- We observe, however, a 15% deviation of the result of this approximation from that of the exact 
calculation for Cs(nw) in the case of the lighter supersymmetric spectrum. In this case, these two results 
seem to have a similar behaviour in 1/ m 2 H ± for any value of to#± , but they are split by terms (m 2 /M| ugY )™, 
resummed in the exact calculation, and, possibly, by intrinsic constants arising from the two-loop calculation. 

What is the impact of this deviation, and of our exact calculation in general, for the BK(B — ► X s j) 
and for an exclusion plot of the charged-Higgs mass remains to be seen. We have not attempted to draw 
such a plot, since a calculation of the W /± contribution at the same precision on m 2 / M| USY is still not yet 
available. Moreover, in the region of supersymmetric parameter space in which the difference between our 
exact calculation and the approximated one of Ref. [1 1 j is largest, one expect also a rather big contribution 
from the chargino-stop exchange. We postpone the presentation of such a plot to later work. 

5 Conclusion 

We have studied the SUSY-QCD corrections to the charged-Higgs contributions to the decay B — * X s "f. 
They are induced by gluino-squark subloops of 0(a s tan/3) and are therefore potentially large in the large- 
tan/3 regime. 

The resulting two-loop diagrams had been dealt in the past in an approximate way. In particular, they 
had been treated in the approximation of an effective Lagrangian with two-Higgs doublets, in which only the 
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nondccoupling operators had been included. For the charged-Higgs contributions to b — > and b — > sg, this 
means that terms of 0(mf, m? H ± /-Mf USY ) or higher, (with MgusY one of the typical squark/gluino masses) 
are neglected. Terms of this type are in general induced by higher-dimensional operators. A truncation of 
the basis of operators in the 2HD effective Lagrangian, is acceptable in particular directions of parameter 
space, in which m H ± <C Msusy- This condition, however, is not generically supported by the most known 
and studied models of supersymmetry breaking. In these, the charged-Higgs mass tends to align along the 
gluino mass, in turn of the same order of the /i parameter. (There are, however, directions in which a special 
tuning among the different supersymmetric parameters allows considerably lower values of m H ±.) 

For the charged-Higgs contributions to b — > sj and b — > sg, we have, therefore included all sublead- 
ing terms, up to ©((m^ , m^ ± /M| USY ) 2 ), obtained through the technique of Heavy Mass Expansion of 
multi-loop integrals. This implicitly takes into account the contribution of all relevant higher-dimensional 
operators that should be added to the 2HD effective Lagrangian. We have also performed the exact calcula- 
tion of all two- loop diagrams correcting at order 0(a s tan/3) the charged-Higgs contribution to b — ► S7 and 
b — > sg. Thus, we have compared to this exact result the different approximations obtained in the HME, i.e. 
1) that in which only the nondecoupling part of the two-loop integrals is retained; 2) that in which terms of 
0(m1 , rn 2 H ± /-Mf USY ) are also added; and finally 3) that in which terms up to 0((toj , m^ ± /M| USY ) 2 ) are 
included. 

We have found that for considerably lighter than the remaining supersymmetric particles, the 
result from the exact calculation and from the three approximations deviate very little one from the other, 
for the coefficient CV(/ziy)- For the coefficient Cs(nw) we find a deviation of 15% between the result of 
approximation 1 ) and that of the exact calculation, if the supersymmetric spectrum is not particularly heavy. 
The calculation presented here is a first step towards the exact evaluation of all supersymmetric contributions 
to B — > X s ^f at order 0(a s tan/3). Only after the completion of such a program, will phenomenological 
analyses be performed and exclusion plots for the masses of the charged-Higgs and other supersymmetric 
particles be attempted. 
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